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Abstract 

In-depth understanding of the retarded oxidation phenomenon observed during the oxidation of 

C^ ', sihcon nanostructures is proposed. The wet thermal oxidation of various sihcon nanostructures 

g 

I ' such as nanobeams, concave/convex nanorings and nanowires exhibits an extremely different and 

C ■ complex behavior. Such effects have been investigated by the modeling of the mechanical stress 

^ generated during the oxidation process explaining the retarded regime. The model describes the 

oxidation kinetics of silicon nanowires down to a few nanometers while predicting reasonable and 

I :«.,,.-.„..,„„...„„.„,,.,-...- ™.. .. .^. .... 

0^ . in silicon oxide through plastic flow. 
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Retarded oxidation where the oxide growth slows down very rapidly with oxidation du- 



liiJs 



ration or with the silicon nanoobject dimension is still a puzzling physical effect [l|-|3|]. This 

physical effect can be viewed as a technological nanoscale tool able to control the nanoob- 

jects shape, size distribution interface properties and could be used in many applications 

J] . However, only very few studies have been dedicated to the understanding of the phe- 



Blii. 



In this work, oxidation kinetics have 



nomenon which remains fragmented and limited 
been investigated both on the experimental and theoretical counterparts in order to improve 
the understanding of the mechanisms of retarded mechanisms and to quantify the amount 
of stress generated at the Si/Si02 interface in silicon nanostructures. 

With the current top-down fabrication capabilities, etched silicon nanostructures includ- 
ing nanobeams, nanorings and nanowires have been fabricated with a high resolution [8| 
and then wet oxidized at 850°C. The seminal work of Kao et al. [9|, [lO| with micrometer 
size of 2D cylindric object is revisited but in the nanometric range. Experimentally, the 
oxidation kinetics have been observed to be strongly dependent on the size and the geom- 
etry of the nanoobject. Fig. [Ha) summarizes the evolution of the oxide thickness as a 
function of the oxidation duration in the case of convex (SiNWs) and concave (Si nanor- 
ings) structures. The oxide growth rate is strongly limited with the oxidation time but is 
faster in a convex structure than in a concave one. The influence of the geometrical effect 
is stronger with smaller inner radius (i.e. 70 nm compared to 430 nm). For the convex 
case, a higher oxide growth rate is related to the larger radii. Then, in order to investigate 
experimentally the influence of silicon nanostructure dimension, nanobeams and nanowires 
of height 240 nm have been oxidized for 10 and 20 minutes. As shown in insert [Tlb).(l), 
a non-uniform oxide growth is classically observed along the sidewall of the beam due to 
the great influence of the top and bottom corners corresponding to a convex and concave 
structure respectively. An oxidized one-dimensional nanostructure with diameters from 40 
nm to 140 nm demonstrated completely different shapes as shown in insert [T]b). (2) with the 
presence of a pinching effect at the bottom of Si nanobeam structures. In order to compare 
the oxidation between the two structures, the oxide thicknesses have been measured in the 
middle of these vertical structures described in SEM images and are plotted in Fig. [Hb). 
The oxide growth on Si nanobeams of previous width L is clearly thicker than a SiNWs 
with the diameter d = L. A size dependent oxidation kinetic was not observed in these 
structures whatever the nanobeam width considered. These experimental results illustrate 



that the sihcon oxidation retarded mechanism is strongly dependent at the nanoscale level 
on the nanoobject i) dimension ii) size and iii) shape. 

These dependences cannot be explained by the standard Deal and Grove oxidation model 



ll| . as for example, a larger oxidant concentration for the smallest particles should in 
principle lead to a higher oxidation rate. Two main theories have been put forward to explain 
the retarded/self-limiting kinetics factor. The first one is the "stress limited reaction rate" 
assumption |6| with a radial stress build-up at the Si/Si02 interface assumed to be linearly 
dependent of the oxidation time up to a critical stress estimated to a few GPa where the 
oxidation rate would be completely negligible. The second theory is the "diffusion limited 
mechanism" associated to a significant increase in the activation energy of the oxidant 
diffusivity in the highly stressed region is put forward |2|. In this framework, the origin of 
self-limited effects would be the oxidant species supply at the interface. Despite the fact that 
an unknown and uncontrolled amount of strain is introduced, no quantitative determination 
of the mechanical stress build-up during processing is provided by the two approaches. 
In order to model the oxidation of cylinder nanoobjects, the extended Deal and Grove 



model in cylindrical coordinates has been used [lO|. The wet oxidation rate v at the Si/Si02 
interface is given by (Eq. [1]) : 

{a-l)C* 1 ^^^ 



N 1 a 6 

with respectively a (b) the inner (outer) radius, a the volume expansion factor of silicon to 
oxide conversion (2.25), A^ the number of oxidant molecules incorporated into a unit volume 
of silicon oxide, the + (resp. -) sign denotes the convex (concave) surface. This equation 
classically takes into account that the surface curvature influences the oxidant concentration 
and in the convex (concave) configuration, the concentration increases (decreases). In our 
approach, both the reaction rate k^^ at the Si/Si02 interface and the diffusivity in the 
silicon oxide Dsi02 ^^^ stress dependent. The reaction rate kg^ is directly proportional to 
the linear rate constant {B/A)[iiq]{T) (5.18 10^°^ nm/s at 850°C) defined in the Deal and 
Grove approach by introducing C* the oxidant solubility in the silicon dioxide and is strongly 
dependent on the radial stress component ar at the Si/Si02 interface : 

fc?. = ^-(5M)[iio](T)exp(^) (2) 



where ks is the Boltzinann constant, T is the oxidation temperature and V^ (15 A^) 
corresponds to the activation volume. A compressive radial stress (o",r < 0) slows down the 
linear rate oxidation rate. The term (S/A)[iio](T) takes into account the influence of the 
[110] crystalline orientation and the factor taking into account orientation effects has been 
calibrated with planar bulk oxidation experiments. Next the oxidant diffusivity, Dsi02 '■ 

Ds.o. = ^-B(T)^^^B„(T).e.,(=^) (3) 

is linearly dependent on the initial parabolic constant Bq{T) (2.68 10^^'^ nm^/s at 850°C) 
and is limited by a compressive (P > 0) hydrostatic pressure P = —0.5 • (a^ + ag) in the 
silicon oxide (Vrf=45 A'^). These assumptions are often estimated to be equivalent to a 



diffusivity dependence with oxide density |l2|. 



A major issue in oxidation modelling is a proper description of the mechanical behavior 



of silicon dioxide as shown, in Fig. |2l and its ability to store or to dissipa 



te mechanical 



energy. A shortcoming is also observed for the viscous standard approach 10 



13| since the 



compressive radial stress at the interface is inversely proportional to the curvature radius of 
the nanoobject and strongly overestimates the stress level [ij]. The main reason is that the 
irreversible atomic rearrangements occurring with large shearing forces [15]] are neglected. 



This plastic flow can been described by a shear dependent viscosity [16[ : 



vir)=Vo{T)-^p—^ (4) 

smh(r/crc) 

where rjoiT) is the low stress viscosity, r is the critical resolved shear stress and cxc is 
the critical stress threshold where plasticity flow should appear (1 GPa). The low stress 
viscosity value (1.4 10^^ Poisse at 850°C) considered is characteristic of a wet oxide with a 



high viscosity induced by the presence of hydroxyl content 17| . Following the expression of 
the critical resolved shear stress T{r) = ^^, it can be underlined that the oxidation growth 
rate (Eq. [1]), the shear dependent viscosity (Eq. [2]) and finally the critical shear stress are 
coupled to each others. The fact that all of these equations must be self-consistently solved 



is often overlooked or not exactly taken into account. Following Rafferty et al. p^, the 
radial (cTj,) and tangential (ag) stress field component in the silicon dioxide of a cylinder 
structure (see Fig. [TJb) can be expressed as : 



aeir) = cr,.(r) - 2r(r) 



(5) 
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with the reduced parameter R = \ — '^ — . Compared to a standard viscous approach 



10| with a constant viscosity, the radial stress build-up has a logarithmic dependence on the 
curvature radius which gives us the opportunity to model the oxidation of cylinder shape 
nanostructures. 

Fig. [3la) shows that the influence of the concave or convex character on the oxidation 
kinetics can be well predicted by the model. As shown in Fig. [3lb) a substantial non-linear 
increase is observed for the compressive radial stress component at the Si/Si02 interface up 
to a few GPa, coupled with an initial tensile hydrostatic pressure (inset Fig. Ob). The radial 
stress build-up is not linear with time as assumed in a previous study [6[. These elements 
clearly indicate that a reaction limited process takes place in the convex configuration. The 
situation is totally different in the concave case where a radial compressive stress build-up 
remains limited whereas the compressive hydrostatic pressure clearly impacts the oxidant 
diffusivity. In that case, the major limiting factor is the diffusion mechanism which reduces 
the oxidant supply. A quasi self-limited oxidation for 70 nm concave structure is observed 
and can be correlated to the occurrence of both a diffusion and reaction limited regime. As 
summarized by table [H the dominant retardation mechanism is strongly dependent on the 
surface shape but can be explained by the variation of the stress field component at the 
Si/Si02 in agreement with previous results [lO|. 

Fig. Hla) presents the linear rate modeling with oxidation duration and SiNW diameters 
reported in Fig[TJb). A strong decrease in the reaction rate with radial stress build-up as 
a function of SiNW diameters is predicted which could be directly correlated to the large 
radial stress build up depicted in Fig. Hlc). Large non-linear build-up as a function of the 
oxidation time and SiNW diameters for the compressive radial stress down to 4 GPa which 
remains compatible with interfacial stress estimated using contact-resonance atomic force 
microscopy [19]. This effect causes the initial retarded effects observed in the oxidation of 
convex nanostructures. These simulations explain well the difference in terms of behavior 
between the nanobeams and the nanowires since finite elements simulations estimate a much 



lower compressive stress build-up in the nanobeam. On the other hand, Fig. |Hb) presents 
the evolution of the parabolic rate with SiNW diameters. A decrease in the parabolic rate is 
observed after a significant delay which can be correlated to the tangential stress relaxation 
as shown in inset of Fig. IHc). This diffusion limited effect generated by a compressive 
hydrostatic pressure is probably much more difficult to control as a time and diameter 
dependence is observed in Fig. Hlb). 

In summary, retarded oxidation kinetics have been investigated at the nanoscale level for 
different silicon nanoobjects. We have demonstrated the influence of the nanoobject dimen- 
sion, size and shape on the oxidation behavior. All of these effects have been correlated to 
the interfacial stress build-up during oxidation. Given the deformation rate level, modelling 
aspects show that plastic relaxation needs to be considered in order to estimate i) a physical 
mechanical stress build-up at the interface and ii) the interface velocity. Both reaction or 
diffusion limited mechanisms must be considered to desccribe retarded oxidation effects in 
silicon nanostructures. 
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TABLE I. Main physical mechanism governing the retarded or self-hmited oxidation with decreas- 
ing convex or concave nanobjects size. 

Character 0"^ P Origin 



Convex /^ <^ Limited reaction rate 
Concave <^ /^ Limited diffusion mechanism 



FIGURES CAPTION 

Fig. [U (a) Convexity vs concavity: the various symbols present the experimental oxide 
thickness for SiNWs (convex nanostructure) with diameters 70 nm, 130 nm and silicon 
nanorings (concave nanostructure) with inner diameters of 70 nm and 430 nm as a function 
of the oxidation time (b) Oxide thickness in the middle edge of the Si nanobeams and SiNWs 
for 10 min and 20 min at 850°C. The experimental trend is described by the dashed lines. 
Inset SEM images: (1) the Si nanobeams (cross-section view), (2) SiNWs (titled view) after 
oxidation where the position set to measure the oxide thickness is reported. 

Fig. 12) Schematics of the concave and convex cylinder nanostructure oxidation and 
resulting stress field in the silicon oxide. The strain in the oxide could be divided into two 
components the deviatoric part associated to shape modification and the dilatational part 
often neglected. Plasticity effects are introduced by considering a non-linear shear dependent 
viscosity //(r). 

Fig. [3] (a) Modelling of the convex/ concave oxidation kinetics using the plastic model, 
b) Evolution of the theoretical radial stress ar at the Si/Si02 interface. Inset provides the 
hydrostatic pressure evolution during the oxidation for the different nanostructures. 

Fig. m (a) Linear rate variation with SiNW diameters showing the impact of the radial 
stress in the initial oxidation regime as the main limiting factor, (b) Highlight of the diffusion 
limited regime which takes place only when a sufficient oxide thickness has been grown, (c) 
Compressive radial stress build-up with decreasing SiNW diameters. Inset: tangential stress 
relaxation. 
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